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A novel information-geometrodynamical approach to chaotic dynamics (IGAC) on curved statis- 
tical manifolds based on Entropic Dynamics (ED) is presented and a new definition of information 
geometrodynamical entropy (IGE) as a measure of chaoticity is proposed. The general classical 
formalism is illustrated in a relatively simple example. It is shown that the hyperbolicity of a non- 
maximally symmetric 6A'^-dimensional statistical manifold A^s underlying an ED Gaussian model 
describing an arbitrary system of SA'^ degrees of freedom leads to linear information-geometric en- 
tropy growth and to exponential divergence of the Jacobi vector field intensity, quantum and classical 
' features of chaos respectively. An information-geometric analogue of the Zurek-Paz quantum chaos 

\ criterion in the classical reversible limit is proposed. This analogy is illustrated applying the IGAC 

, to a set of n-uncoupled three-dimensional anisotropic inverted harmonic oscillators characterized 

^ ' by a Ohmic distributed frequency spectrum. 

I. INTRODUCTION 

, The lack of a unified characterization of chaos in classical and quantum dynamics is well-known. In the Riemannian 
[l| and Finslerian [3| (a Finsler metric is obtained from a Riemannian metric by relaxing the requirement that the 
metric be quadratic on each tangent space) geometrodynamical approach to chaos in classical Hamiltonian systems, 
an active field of research concerns the possibility of finding a rigorous relation among the sectional curvature, the 
' Lyapunov exponents, and the Kolmogorov-Sinai dynamical entropy (i. e. the sum of positive Lyapunov exponents) 
. [3|. The largest Lyapunov exponent characterizes the degree of chaoticity of a dynamical system and, if positive, it 
■ measures the mean instability rate of nearby trajectories averaged along a sufhciently long reference trajectory. More- 
\^ , over, it is known that classical chaotic systems are distinguished by their exponential sensitivity to initial conditions 
t;^ ■ and that the absence of this property in quantum systems has lead to a number of different criteria being proposed 
' for quantum chaos. Exponential decay of fidelity, hypersensitivity to perturbation, and the Zurek-Paz quantum chaos 
criterion of linear von Neumann's entropy growth Q are some examples Q. These criteria accurately predict chaos 
' in the classical limit, but it is not clear that they behave the same far from the classical realm. 

. The present work makes use of the so-called Entropic Dynamics (ED) [g^. ED is a theoretical framework that arises 
' from the combination of inductive inference (Maximum relative Entropy Methods, |J3]) and Information Geometry 
(Riemannian geometry applied to probability theory) (IG) As such, ED is constructed on statistical manifolds. 
It is developed to investigate the possibility that laws of physics - either classical or quantum - might reflect laws of 
^ , inference rather than laws of nature. 

This article is a follow up of a series of the authors works [1, [T^iIHI- In this paper, the ED theoretical framework is 
used to explore the possibility of constructing a unified characterization of classical and quantum chaos. We investigate 
a system with 3iV degrees of freedom (microstates) , each one described by two pieces of relevant information, its mean 
expected value and its variance (Gaussian statistical macrostates) . This leads to consider an ED model on a non- 
maximally symmetric GA^-dimensional statistical manifold AAs- It is shown that M.s possesses a constant negative 
Ricci curvature that is proportional to the number of degrees of freedom of the system, Rm, — SN. It is shown that 
the system explores statistical volume elements on Ms at an exponential rate. We define a dynamical information- 
geometric entropy of the system and we show it increases linearly in time (statistical evolution parameter) and 
is moreover, proportional to the number of degrees of freedom of the system. The geodesies on Aig are hyperbolic 
trajectories. Using the Jacobi-Levi-Civita (JLC) equation for geodesic spread, it is shown that the Jacobi vector 
field intensity Jm, diverges exponentially and is proportional to the number of degrees of freedom of the system. 
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Thus, j and Jms are proportional to the number of Gaussian-distributed microstates of the system. This 
proportionahty leads to conclude there is a substantial link among these information-geometric indicators of chaoticity. 

Finally, an information-geometric analog of the Zurek-Paz quantum chaos criterion is suggested. We illustrate this 
point by use of an n-set of inverted harmonic oscillators (IHO). In the ED formalism, the IHO system is described 
by a curved n-dimensional statistical manifold that is conformally related to an Euclidean one. 



II. SPECIFICATION OF THE GAUSSIAN ED-MODEL 

Maximum relative Entropy (ME) methods are used to construct an ED model that follows from an assumption about 
what information is relevant to predict the evolution of the system. Given a known initial macrostate (probability 
distribution) and that the system evolves to a final known macrostate, the possible trajectories of the system are 
examined. A notion of distance between two probability distributions is provided by IG. As shown in fisj this 
distance is quantified by the Fisher-Rao information metric tensor. 

We consider an ED model whose microstates span a 3iV-dimensional space labelled by the variables |-^| = 
x^^\...., with EE (x[°'\ xi"\ 4"^), a = 1,...., N and xi"^ e R with a = 1, 2, 3. We assume the 

only testable information pertaining to the quantities xi"^ consists of the expectation values (xi"^ ) and variance 



Axi"'' = y y^xi"'' — (^xi"^^^ y. The set of these expectation values define the 6A^-dimensional space of macrostates 

of the system. A measure of distinguishability among the states of the ED model is obtained by assigning a probability 

distribution P (^X to each macrostate 6 where |e| = ^^^'^0["\ ^'^^0["^^ with a = 1, 2,...., TV and a = 1, 2, 3. 

The process of assigning a probability distribution to each state endows A4s with a metric structure. Specifically, the 
Fisher-Rao information metric defined in (jl2p is a measure of distinguishability among macrostates. It assigns an IG 
to the space of states. 



A. The Gaussian statistical manifold Ms 



We consider an arbitrary system evolving over a 3A''-dimensional space. The variables |-^| = x*-^-*,...., x*^^)} 

label the 3iV-dimensional space of microstates of the system. All information relevant to the dynamical evolution 
of the system is assumed to be contained in the probability distributions. For this reason, no other information is 
required. Each macrostate may be viewed as a point of a 6iV-dimensional statistical manifold with coordinates given 



by the numerical values of the expectations '-'^^6'*,"'' = (|xi"''^ and — Axi"'' 

available information is contained in the following 6A^ information constraint equations. 



The 



4-c)o 



dx 



(1) 



Axi") 



dx, 



where ^^^dl"^ = v"/ ^^'^ " ^^"^'^ ^^^^ a = 1, 2,...., iV and a = 1, 2, 3. The probability distributions Pi"^ 

are constrained by the conditions of normalization. 



+ 00 



dxi")pi") (xM|we\(^^eO = i. 



(2) 



The Gaussian distribution is identified by information theory as the maximum entropy distribution if only the expec- 
tation value and the variance are known. ME methods allows to associate a probability distribution P (^X 9 ^ to each 
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point in the space of states 8. The distribution that best reflects the information contained in the prior distribution 
m updated by the information ^^a;i"''^ , Axi"''^ is obtained by maximizing the relative entropy 



e log 



{^} 




(3) 



As a working hypothesis, the prior m is set to be uniform since we assume the lack of prior available information 
about the system (postulate of equal a priori probabilities). Upon maximizing ([3]), given the constraints ([1]) and 



we obtain 



N 3 



^(^|0) = nn^i"^ (4"' 



(4) 



where 



exp 



pa 



(5) 



and (1)61^"^ = Ma"\ (^^^^"^ — ai"'. For the rest of the paper, unless stated otherwise, the statistical manifold Ms will 
be defined by the following expression. 



Ms=[p{x e) in Q : X e M^^, 9 £ = (-oo, + oo)^ x (0, + oo) 



37V 



(6) 



The probability distribution @ encodes the available information concerning the system. Note we assumed uncoupled 
constraints among microvariables xi"'. In other words, we assumed that information about correlations between the 
microvariables need not to be tracked. This assumption leads to the simplified product rule (jj]). However, coupled 
constraints would lead to a generalized product rule in @ and to a metric tensor (fT^ with non-trivial off-diagonal 
elements (covariance terms). For instance, the total probability distribution P (x, y\^J^x^ <^x, l^y, <^y) of two dependent 
Gaussian distributed microvariables x and y reads 



P {x, y\n^, a^, fly, ay) 



27raxO'yVl — r^ 



(7) 



X exp ■ 



1 



2(1-7-2) 

where r G ( — 1, + 1) is the correlation coefficient given by 

((a; - (x)) {y - (y))) (xy) - (x) (y) 



^{x~{x))^{y-{y)) 
The metric induced by ([7]) is obtained by use of p^ . the result being 



(8) 



9^3 





r 

o"xO"!/(f^-l) 


















crxO-„(r2-l) 





(9) 



where i, j = 1, 2, 3, 4. The Ricci curvature scalar associated with manifold characterized by 

8 (r2 -2)+ 2r2 (Sr^ - 2) 



is given by 



R 



— R - — 

1 to 



;(r2 - 1) 



(10) 
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It is clear that in the limit r ^ 0, the ofF-diagonal elements of gij vanish and the Scalar R reduces to the result obtained 
in [Toj . namely R = —2 < 0. Correlation terms may be fictitious. They may arise for instance from coordinate 
transformations. On the other hand, correlations may arise from external fields in which the system is immersed. In 
such situations, correlations among xi"'' effectively describe interaction between the microvariables and the external 
fields. Such generalizations would require more delicate analysis. Before proceeding, a comment is in order. Most 
probability distributions arise from the maximum entropy formalism as a result of simple statements concerning 
averages (Gaussians, exponential, binomial, etc.). Not all distribution are generated in this manner however. Some 
distributions are generated by combining the results of simple cases (multinomial from a binomial) while others 
are found as a result of a change of variables (Cauchy distribution). For instance, the WeibuU and Wigner-Dyson 
distributions can be obtained from an exponential distribution as a result of a power law transformation [l5| . 



Metric structure of Ms 



We cannot determine the evolution of microstates of the system since the available information is insufficient. Not 
only is the information available insufficient but we also do not know the equation of motion. In fact there is no 
standard "equation of motion". Instead we can ask: how close are the two total distributions with parameters {t^a^\ 
a'^a'^) and (/ii""* + d^'"a'\ ui"'' + da\^^)l Once the states of the system have been defined, the next step concerns the 
problem of quantifying the notion of change from the state O to the state 8 + dQ. A convenient measure of change 



is distance. The measure we seek is given by the dimensionless distance ds between P (^X 8^ 

ds"^ = g^^.dQ^'dQ" with ii,v=l, 2,.., 6A^, 



and P X 



Q + dQ 



where 



g^,. ^ / dXPiX 



8 



d\ogP[X 



8 dlogPiX 
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(11) 



(12) 



is the Fisher-Rao information metric. Substituting @ into (|12p . the metric g^i, on M.s becomes a 6A^ x 6A^ matrix 
M made up of 3 A'' blocks M2X2 with dimension 2x2 given by. 



M2 







2 X (ai")) 

with a = 1, 2,...., A^ and a = 1,2,3. From the "length" element ((TTl) reads. 



(13) 



N 3 



a — 1 a— 1 



1 



Ac) 



Ml 



(q)2 



(14) 



We bring attention to the fact that the metric structure of Ais is an emergent (not fundamental) structure. It arises 
only after assigning a probability distribution P [X 8 ) to each state 8. 



2. Curvature of Ms 

Given the Fisher-Rao information metric, we use standard differential geometry methods applied to the space of 
probability distributions to characterize the geometric properties oi Ms- Recall that the Ricci scalar curvature R is 
given by, 

R^g^'^'R^,, (15) 



where g^^gvp — 6^ so that g^" — {gp^u) ^- The Ricci tensor is given by. 



Rpu — d^^'li, d^T^x + T^nu^l„ '^Zl^'i^ir (16) 
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The Christoffel symbols appearing in the Ricci tensor are defined in the standard manner as, 

= (9^5.. + d.9i.a - d,g^,) . (17) 
Using (|13p and the definitions given above, we can show that the Ricci scalar curvature becomes 

The scalar curvature is the sum of all sectional curvatures of planes spanned by pairs of orthonormal basis elements 
{ep = <9ep(p)} of the tangent space TpMs with p S Ms, 

K (a, b) = R^'^P-'^'^b'^'^'b- , a = y (a, h") e^, (19) 

i9^.ag.p-g^.p9.a)aH-aPb-' ^ 

where (cp, h'^) — (5^. Notice that the sectional curvatures completely determine the curvature tensor. From (fTSl) we 
conclude that Ais is a GA^-dimensional statistical manifold of constant negative Ricci scalar curvature. A detailed 
analysis on the calculation of Christoffel connection coefficients using the ED formalism for a four-dimensional manifold 
of Gaussians can be found in [l3| . 



3. Anisotropy and Compactness 

It can be shown that is not a pseudosphere (maximally symmetric manifold). The first way this can be 
understood is from the fact that the Weyl Projective curvature tensor fT6l| (or the anisotropy tensor) W^^pa defined 

by 

Wp^^pcr = Rfj,vpa ~ ^ ^ (Sfcr^MP 9vp9li<y) i (20) 

with n ~ 6N in the present case, is non- vanishing. In (j20p . the quantity R^ypa is the Riemann curvature tensor 
defined in the usual manner by 

Ppty = d^r^pp ~ dpT^p^ + r" \cr^^ pp — r" xpV^ p^. (21) 

Considerations regarding the negativity of the Ricci curvature as a strong criterion of dynamical instability and the 
necessity of compactness of in "true" chaotic dynamical systems is under investigation (l3 |. 

The issue of symmetry of A^s can alternatively be understood from consideration of the sectional curvature. In 
view of pop , the negativity of the Ricci scalar implies the existence of expanding directions in the configuration space 
manifold M.s- Indeed, from (|18p one may conclude that negative principal curvatures (extrema of sectional curvatures) 
dominate over positive ones. Thus, the negativity of the Ricci scalar is only a sufficient (not necessary) condition 
for local instability of geodesic fiow. For this reason, the negativity of the scalar provides a strong criterion of local 
instability. Scenarios may arise where negative sectional curvatures are present, but the positive ones could prevail in 
the sum so that the Ricci scalar is non-negative despite the instability in the fiow in those directions. Consequently, 
the signs of the sectional curvatures are of primary significance for the proper characterization of chaos. 

Yet another useful way to understand the anisotropy of the AAg is the following. It is known that in n dimensions, 
there are at most "■^^^'^'^ independent Killing vectors (directions of symmetry of the manifold). Since A^s is not a 
pseudosphere, the information metric tensor does not admit the maximum number of Killing vectors Ki, defined as 

Ck9^u = Dp,K, + D,Kp = 0, (22) 

where D^, given by 

DpK, - d^,K, - TP^Kp, (23) 



is the covariant derivative operator with respect to the connection F defined in (jl7p . The Lie derivative CKg^u of the 
tensor field gp^i, along a given direction K measures the intrinsic variation of the field along that direction (that is, 
the metric tensor is Lie transported along the Killing vector) [iTj . Locally, a maximally symmetric space of Euclidean 
signature is either a plane, a sphere, or a hyperboloid, depending on the sign of R. In our case, none of these scenarios 
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occur. As will be seen in what follows, this fact has a significant impact on the integration of the geodesic deviation 
equation on Mg. At this juncture, we emphasize it is known that the anisotropy of the manifold underlying system 
dynamics plays a crucial role in the mechanism of instability. In particular, fluctuating sectional curvatures require 
also that the manifold be anisotropic. However, the connection between curvature variations along geodesies and 
anisotropy is far from clear and is currently under investigation. 

Krylov was the first to emphasize the use of i? < as an instability criterion in the context of an iV-body 
system (a gas) interacting via Van der Waals forces, with the ultimate hope to understand the relaxation process in a 
gas. However, Krylov neglected the problem of compactness of the configuration space manifold which is important 
for making inferences about exponential mixing of geodesic flows 19]. Why is compactness so significant in the 
characterization of chaos? True chaos should be identified by the occurrence of two crucial features: 1) strong 
dependence on initial conditions and exponential divergence of the Jacobi vector field intensity, i.e., stretching of 
dynamical trajectories; 2) compactness of the configuration space manifold, i.e., folding of dynamical trajectories. 
Compactness 0, H^l is required in order to discard trivial exponential growths due to the unboundedness of the 
"volume" available to the dynamical system. In other words, the folding is necessary to have a dynamics actually 
able to mix the trajectories, making practically impossible, after a finite interval of time, to discriminate between 
trajectories which were very nearby each other at the initial time. When the space is not compact, even in presence 
of strong dependence on initial conditions, it could be possible in some instances (though not always), to distinguish 
among different trajectories originating within a small distance and then evolved subject to exponential instability. 

The statistical manifold defined in © is compact. This can be seen as follows. It is known form IG that there 
is a one-to-one relation between elements of the statistical manifold and the parameter space. More precisely, the 
statistical manifold is homeomorphic to the parameter space Vq. This implies the existence of a continuous, 
bijective map /iaIs, ©ei 



Ms^ P(x ^Ve (24) 



where hjj^ = P (^X O ^ . The inverse image hjj^ is the so-called homeomorphism map. In addi- 

tion, since homeomorphisms preserve compactness, it is sufficient to restrict ourselves to a compact subspace of the 
parameter space Vq in order to ensure that A^s is itself compact. 

III. CANONICAL FORMALISM FOR THE GAUSSIAN ED-MODEL 

The geometrization of a Hamiltonian system by transforming it to a geodesic flow is a well-known technique of 
classical mechanics associated with the name of Jacobi [21j. Transformation to geodesic motion is obtained in two 
steps: 1) conformal transformation of the metric; 2) rescaling of the time parameter [22|. The reformulation of 
dynamics in terms of a geodesic problem allows the application of a wide range of well-known geometrical techniques 
in the investigation of the solution space and properties of equations of motions. The power of the Jacobi reformulation 
is that all of the dynamical information is collected into a single geometric object - the manifold on which geodesic 
flow is induced - in which all the available manifest symmetries are retained. For instance, integrability of the system 
is connected with the existence of Killing vectors and tensors on this manifold [1^ [13] ■ 

In this section we study the trajectories of the system on Aig- We eniphasize ED can be derived from a standard 
principle of least action (of Maupertuis-Euler-Lagrange- Jacobi type) [y, |2^. The main differences are that the 
dynamics being considered here, namely ED, is defined on a space of probability distributions A4s, not on an ordinary 
linear space V and the standard coordinates of the system are replaced by statistical macrovariables 8^. The 
geodesic equations for the macrovariables of the Gaussian ED model are given by, 

^ + C,^^^ = (25) 

with /i = 1, 2,..., 6A^. Observe the geodesic equations are nonlinear second order coupled ordinary differential equa- 
tions. They describe a reversible dynamics whose solution is the trajectory between an initial and a final macrostate. 
The trajectory can be equally well traversed in both directions. 

A. Geodesies on Ms 

We determine the explicit form of for the pairs of statistical coordinates {^J.^"\ (ri"-*). Substituting the expres- 
sion of the Christoffel connection coefficients into (|25p. the geodesic equations for the macrovariables ^i"'* and cri"'* 
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associated to the microstate cci""* become, 

with a = 1, 2,...., A'^ and a = 1, 2, 3. This is a set of coupled ordinary differential equations, whose solutions are 



(27) 



cosh(2/3<"'r)-sinh(2/3<°'r)+i^-^ 



J") _ cosh(/j(°)r)-sinh(/3i°)r) 

cosh (2/3i°' t) -sinh (2/3i°' r ) + ^ ° / .j 

8(/3a ) 

The quantities Bi"\ ci"\ (3^^^ are real integration constants that can be evaluated upon specification of boundary 
conditions. We are interested in the stability of the trajectories on A^g. It is known f25'| that the Riemannian 
curvature of a manifold is intimately related to the behavior of geodesies on it. If the Riemannian curvature of a 
manifold is negative, geodesies (initially parallel) rapidly diverge from one another. For the sake of simplicity, we 

assume very special initial conditions: = A, (3^"^ = A S K"*", Cq""* = 0, Va = 1, 2,...., TV and a = 1, 2, 3. However, 
the conclusions drawn can be generalized to more arbitrary initial conditions. We observe that since every maximal 
geodesic is well-defined for all temporal parameters t, M.s constitute a geodesically complete manifold [1^. It is 
therefore a natural setting within which one may consider global questions and search for a weak criterion of chaos 
i. 

IV. EXPONENTIAL DIVERGENCE OF THE JACOBI VECTOR FIELD INTENSITY 

The actual interest of the Riemannian formulation of the dynamics stems form the possibility of studying the 
instability of natural motions through the instability of geodesies of a suitable manifold, a circumstance that has 
several advantages. First of all a powerful mathematical tool exists to investigate the stability or instability of a 
geodesic flow: the Jacobi-Levi-Civita equation for geodesic spread [27,]. The JLC-equation describes covariantly how 
nearby geodesies locally scatter. It is a familiar object both in Riemannian geometry and theoretical physics (it 
is of fundamental interest in experimental General Relativity). Moreover the JLC-equation relates the stability or 
instability of a geodesic flow with curvature properties of the ambient manifold, thus opening a wide and largely 
unexplored held of investigation of the connections among geometry, topology and geodesic instability, hence chaos. 

Consider the behavior of the one-parameter family of neighboring geodesies ^Gms ('^) = {®!t( ('''! ^)}AeR+' 
where 

mI"^ (r; A) 



A2 1 



2A eosh (2At) - sinh (2At) - 

^(a) ^ cosh (At) - sinh (At) 

" cosh (2At) -sinh (2At) + 



A2 

8A^ 



A2 



(28) 



with a = 1, 2,...., N and a = 1, 2, 3. The relative geodesic spread on a (non-maximally symmetric) curved manifold 
as A^s is characterized by the Jaeobi-Levi-Civita equation, the natural tool to tackle dynamical chaos [itI. [27|. 

_ + fl.^_,e»^=0 ,29) 

where the Jacobi vector held J** is deflned as. 



V d\ 



SX. (30) 

r— const 



Notice that the JLC-equation appears intractable already at rather small N. For isotropic manifolds, the JLC-equation 
can be reduced to the simple form, 

^^+KJ'^ = 0, (,= !,...., 6N (31) 
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where K is the constant value assumed throughout the manifold by the sectional curvature. The sectional curvature 
of manifold A^^ is the 6A^-dimensional generalization of the Gaussian curvature of two-dimensional surfaces of R'^. If 
K <0, unstable solutions of equation ([3T|l assumes the form 



J(r) 



1 



=uj (0) sinh I 



(32) 



once the initial conditions are assig ned as J (0) = 0, = uj (0) and K < 0. Equation ^ forms a system of 6A^ 

coupled ordinary differential equations linear in the components of the deviation vector field (|30p but nonlinear in 
derivatives of the metric ([H]). It describes the linearized geodesic flow: the linearization ignores the relative velocity of 
the geodesies. When the geodesies are neighboring but their relative velocity is arbitrary, the corresponding geodesic 
deviation equation is the so-called generalized Jacobi equation 28, '29] . The nonlinearity is due to the existence of 
velocity-dependent terms in the system. Neighboring geodesies accelerate relative to each other with a rate directly 
measured by the curvature tensor Rafj-yS- Substituting (j28p in (|29p and neglecting the exponentially decaying terms 
in 50'^ and its derivatives, integration of leads to the following asymptotic expression of the Jacobi vector field 
intensity, 



J 



Ms 



\\J\\ = {g^^.J^ry- 



3Ne 



At 



(33) 



We conclude that the geodesic spread on A^^ is described by means of an exponentially divergent Jacobi vector field 
intensity Jm,, a classical feature of chaos. In our approach the quantity Aj, 



Aj =^ lim — In 

T — *oo T 



(34) 



would play the role of the conventional Lyapunov exponents. 



V. LINEARITY OF THE INFORMATION GEOMETRODYNAMICAL ENTROPY 

We investigate the stability of the trajectories of the ED model considered on Ms- It is known [IH] that the 
Riemannian curvature of a manifold is closely connected with the behavior of the geodesies on it. If the Riemannian 
curvature of a manifold is negative, geodesies (initially parallel) rapidly diverge from one another. For the sake of 

simplicity, we assume very special initial conditions: Bi"^ = A, f]'-"^ = A € R+, ci"^ = 0, V a = 1, 2,...., N and 
a = 1, 2, 3 . However, the conclusion we reach can be generalized to more arbitrary initial conditions. Recall Ms 

is the space of probability distributions |p (^X | labeled by 6A^ statistical parameters Q. These parameters are 

the coordinates for the point P, and in these coordinates a volume element dVM^ reads, 

N 3 /2 

dVMs - ^/5rf'^© = n n -T^h^dfii-Uai^^ . (35) 
The volume of an extended region AV^^ (r; A) of Aig is defined by, 

A^M. (r; A) n n J^d^i^'dai-^ (36) 

where /li"-* (r) and ai"'' (r) are given in (|28p . The quantity that encodes relevant information about the stability of 
neighboring volume elements is the the average volume Vm, (t; A), 

T 

Vm. (r; A) ee {AVm. (t; A))|^ ^ J AVm. (t'; A) dr' g^a^At^ (37) 





This asymptotic regime of diffusive evolution in (I37|) describes the exponential increase of average volume elements on 
Ms- The exponential instability characteristic of chaos forces the system to rapidly explore large areas (volumes) of the 
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statistical manifold. It is interesting to note that this asymptotic behavior appears also in the conventional description 
of quantum chaos where the entropy increases linearly at a rate determined by the Lyapunov exponents ^30i] . The 
linear increase of entropy as a quantum chaos criterion was introduced by Zurek and Paz In our information- 
geometric approach a relevant quantity that can be useful to study the degree of instability characterizing the ED 
model is the information-geometric entropy defined as, 

Sm/= li^n logVMjT; \) . (38) 

r — >oo 

Substituting ^ in (|38p . we obtain 
Sm^ = lim log 

r — ^oo 

The entropy Sm, in ([39|) is the asymptotic limit of the natural logarithm of the statistical weight (AV^i^),- defined 
on Ms- Its linear growth in time is reminiscent of the aforementioned quantum chaos criterion. Indeed, equation (|39p 
may be considered the information-geometric analog of the Zurck-Paz chaos criterion. 
In conclusion, we have shown, 

Rm. = -3iV, Sm. « SiVAr, Jms ~ 3Ne^\ (40) 

The Ricci scalar curvature Rm^i the information-geometric entropy ^^^d the Jacobi vector field intensity Jms 
are proportional to the number of Gaussian-distributed microstates of the system. This proportionality leads to the 
conclusion that there exists a substantial link among these information-geometric measures of chaoticity, namely 

RMs ^ ^ Jms- (41) 

Equation (|4ip , together with the information-geometric analog of the Zurek-Paz quantum chaos criterion, equation 
((39)) . represent the fundamental results of this work. We believe our theoretical modelling scheme may be used to 
describe actual systems where transitions from quantum to classical chaos scenario occur, but this requires additional 
analysis. In the following section, we briefly consider some similarities among the von Neumann, Kolmogorov-Sinai 
and Information-Geometrodynamical entropies. 



N 



nn 



a=la=l (0) 



'(0) 



dr' 



3NX7 



(39) 



VI. ON THE VON NEUMANN, KOLMOGOROV-SINAI AND INFORMATION 
GEOMETRODYNAMICAL ENTROPIES 



In conventional approaches to chaos, the notion of entropy is introduced, in both classical and quantum physics, as 
the missing information about the systems fine-grained state [3l| . For a classical system, suppose that the phase 
space is partitioned into very fine-grained cells of uniform volume Av, labelled by an index j. If one does not know 
which cell the system occupies, one assigns probabilities pj to the various cells; equivalently, in the limit of infinitesimal 
cells, one can use a phase-space density p (Xj) = Then, in a classical chaotic evolution, the asymptotic expression 
of the information needed to characterize a particular coarse-grained trajectory out to time r is given by the Shannon 
information entropy (measured in bits), 

^ciasticii^ =- fdXp (X) log, ip (X) Av)=-Y^ p, log, p, ^ K.T. (42) 

where p{X) is the phase-space density and pj = -g;^ is the probability for the corresponding coarse-grained trajectory. 

•^ciassica,? ^hc missiug information about which fine-grained cell the system occupies. The quantity K. represents 
the linear rate of information increase and it is called the Kolmogorov-Sinai entropy (or metric entropy) (/C is the 
sum of positive Lyapunov exponents, /C = "^j^j )■ K. quantifies the degree of classical chaos. It is worthwhile 
emphasizing that the quantity that grows asymptotically as /Cr is really the average of the information on the left side 
of equation (|^^ . This distinction can be ignored however, if we assume that the chaotic system has roughly constant 
Lyapunov exponents over the accessible region of phase space. In quantum mechanics the fine-grained alternatives 
are normalized state vectors in Hilbert space. From a set of probabilities for various state vectors, one can construct 
a density operator 

P^T.^^ I^j) ' P = l^j) • (43) 
j 
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The normalization of the density operator, tr (p) = 1, imphes that the eigenvalues make up a normalized probability 
distribution. The von Neumann entropy of the density operator p (measured in bits) |32i | , 

'5'qutntum =-tr{p l0g2 P) = " l^gs (44) 

j 

can be thought of as the missing information about which eigenvector the system is in. Entropy quantifies the degree 
of unpredictability about the system's fine-grained state. 

Recall that decoherence is the loss of phase coherence between the set of preferred quantum states in the Hilbert 
space of the system due to the interaction with the environment. Moreover, decoherence induces transitions from 
quantum to classical systems. Therefore, classicality is an emergent property of an open quantum system. Motivated 
by such considerations, Zurek and Paz investigated implications of the process of decoherence for quantum chaos. 

They considered a chaotic system, a single unstable harmonic oscillator characterized by a potential V (x) = — 
(A is the Lyapunov exponent), coupled to an external environment. In the reversible classical limit |33| |. the von 
Neumann entropy of such a system increases linearly at a rate determined by the Lyapunov exponent, 

Cantu::^ (Zurek-Paz) Ar. (45) 

Notice that the consideration of 3N uncoupled identical unstable harmonic oscillators characterized by potentials 
Vi (x) = ^ (Ai = Aj-; i, j = 1, 2,..., 3A^) would simply lead to 

S^;^ (Zurek-Paz) 3iVAr. (46) 

The resemblance of equations ([5^ and (|46p is remarkable. In what follows, we apply our information geometrical 
method to an n-set (for n — 2) of uncoupled inverted harmonic oscillators, each with different frequency, and show we 
obtain asymptotic linear IGE growth. The case for arbitrary n-set in three dimensions is presented in the Appendix. 



VII. THE INFORMATION GEOMETRY OF A 2-SET OF INVERTED HARMONIC OSCILLATORS 

(IHO) 

In this section, our objective is to characterize chaotic properties of a 2-set of one-dimensional inverted harmonic 
oscillators, each with different frequency lui ^ uj2 using the formalism presented in this paper. We will study the 
asymptotic behavior of the geometrodynamical entropy and the functional dependence of the Ricci scalar curvature 

(2) 

of the 2-dimensional manifold M j/iq underlying the ED model of the IHOs on the frequencies , i = 1, 2. Recent 
investigation explore the possibility of using well established principles of inference to derive Newtonian dynamics from 
relevant prior information codified into an appropriate statistical manifold [s^l- In that work the basic assumption 
is that there is an irreducible uncertainty in the location of particles so that the state of a particle is defined by a 
probability distribution. The corresponding configuration space is a statistical manifold the geometry of which is 
defined by the information metric. The trajectory follows from a principle of inference, the method of Maximum 
Entropy. There is no need for additional "physical" postulates such as an action principle or equation of motion, 
nor for the concept of mass, momentum and of phase space, not even the notion of time. The resulting "entropic" 
dynamics reproduces Newton's mechanics for any number of particles interacting among themselves and with external 
fields. Both the mass of the particles and their interactions are explained as a consequence of the underlying statistical 
manifold. 

In what follows, we introduce the basics of the general formalism for an n-set of IHOs. This approach is sim- 
ilar (mathematically but not conceptually) to the geometrization of Newtonian dynamics used in the Riemannian 
geometrodynamical to chaos [l|, [sH] ■ 



A. Informational geometrization of Newtonian dynamics 

Newtonian dynamics can be recast in the language of Riemannian geometry applied to probability theory, namely. 
Information Geometry. In our case, the system under investigation has n degrees of freedom and a point on the 
n dimensional configuration space manifold -M^PjIq is parametrized by the n Lagrangian coordinates (^i,...., (?„). 
Moreover, the system under investigation is described by the Lagrangian £, 

c^T U,..Aa) - <f {ei,.; 0„) = ls,,0,e, + IJ2^P^j (47) 

V / j=i 
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so that the Hamiltonian function H. =T + <I> = i? is a constant of motion. For the sake of simpUcity, let us set i? = 1. 
According to the principle of stationary action - in the form of Maupertuis - among all the possible isoenergetic paths 
7 [t) with fixed end points, the paths that make vanish the first variation of the action functional 

dC ■ 

— e^dr (48) 

are natural motions. As the kinetic energy T is a homogeneous function of degree two, we have 2T ~ Oi — , and 

ae, 

Maupertuis' principle reads 



/ 2TdT = 0. (49) 



51 = 5 

The manifold Mj^q is naturally given a proper Riemannian structure. In fact, let us consider the matrix 

(7„ (0i,.., ^?„) = [l-$(0i,.., 0„)]J,y (50) 

so that Maupertuis' principle becomes 

5 f Tdt = 5 f {T^)^dT = 5 f \[l-^{0i,..,9^)]5,,e,ejy 

= 5 I ig^j0^9■] \t = 5 [ ds = Q,ds^ = g^jdO'dO^ (51) 

thus natural motions are geodesies of M^f^o^ pi'ovided we define ds as its arclength. The metric tensor gj (•, •) of 
■^/ffo ^^^^ defined by 

g ^ g,jde^ ® dO^ (52) 

where (d6'\ , ^6*") is a natural base of TqM^j^q - the cotangent space at the point 6* - in the local chart (6*^,...., 0") . 

This is known as the Jacobi metric (or kinetic energy metric). Denoting by V the canonical Levi-Civita connection, 
the geodesic equation 
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becomes, in the local chart (9^ 



V^7 = (53) 



, , 



ds^ ds ds 



where the Christoffel coefficients are the components of V defined by 



rj, = {d0\ V,ek) = iff™ {d.gum + ^kg„^J - d^g^u) , (55) 



with di = Since gtj {9i,--, On) = [1 — $ (6*1,.., On)] 5ij, from the geodesic equation we obtain 



d(l-^)dO^ dO' ,,, a ( 1 - $) dO'' dO 
2 ^ J — 1 9 9k 



dOj ds ds dOj ds ds 



ds^ 2(1-$ 

whereupon using ds^ = (1 — $)^ dr^, we verify that ([5^)1 reduces to 

d^o' , a$(0i,.., 0„) 



(56) 



dr^ dO^ 



0, i = l,..., n. (57) 



Equation (I57p are Newton's equations. It is worthwhile emphasizing that the transformation to geodesic motion on 
a curved statistical manifold is obtained in two key steps: the conformal transformation of the metric^ 5ij — s- gij = 
(1 — $) 5ij and, the rescaling of the temporal evolution parameter, dr"^ —> ds"^ = 2 (1 — $)^ dr^ . 
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B. The 2-set of inverted anisotropic one-dimensional harmonic oscillators 



As a simple physical example, we examine the IG associated with a 2-set configuration of IHOs. In this case, the 
metric tensor g^j appearing in (jSOp takes the form 



g^3 [Ol, 02) = [1 - * (01, 02)] ■ % {9i, 02) with Z, J = 1, 2. 

where the function $ 02) is given by, 

2 



<^{9i,02) = Y.'^j{0,),<^j{0,) 



(58) 



(59) 



Hence the metric tensor gij on M j^q becomes, 



i + i(^?^? + ^i^2) 



,2/)2 



,2/)2\ 



Using the standard definition of the Ricci scalar ((T5)) . we obtain 

4(0?C^l+(?^^|)-4((??+0^)^2^2_8(^2+^2) 
-'l,.(2) (ti'l, Wo) = 5 . 

{0livl + 0livl + 2)' 

In the limit of a flat frequency spectrum, cji = 0^2 = tj, the scalar curvature (|6ip is constantly negative, 

-16^2 



_ {^) 



[2 + (0? + 0l) 



< 0, Vcj > 0. 



(60) 



(61) 



(62) 



However, in presence of distinct frequency values, wi 7^ lu2, it is possible to properly choose the oj's so that 
-R. .(2) {loi, UJ2) becomes either negative or positive. In addition, we notice that the manifold underlying the IHO 
model is anisotropic since its associated Weyl projective curvature tensor components are non-vanishing. For the 
special case, toi = 0J2, we obtain 



W1212 (t^) 



8cj4 {el + el) + 2uj^ {e\ + el) + Auj^e 



1^2 



=)2,,,2 I fl2,^,2 



(63) 



(2) 

Clearly, the frequency parameter u) drives the degree of anisotropy of the statistical manifold M j^o ^ expected, 
in the limit of vanishing w, we recover the flat (i? = 0), isotropic {W = 0) Euclidean manifold characterized by metric 
5ij. This result is a concrete example of the fact that conformal transformations change the degree of anisotropy 
of the ambient statistical manifold underlying the Newtonian dynamics. Our only remaining task is to compute the 
information geometrodynamical entropy S ..(2) (r; Wi, ^2)1 defined as 



Sm(^) {t; wi, W2) =^ lim log (AT/ (2) (r; wi, ^2)) • 
The quantity ( AV". .(2) (t; lji, UJ2)) appearing in ([64ll is the average volume element, defined by 

T 

(Ay (2) (T;a;i,W2)) ^ - AV (2) (t; ui, UJ2) dr' , 



with the statistical volume element AV..{2) given by 

■'^ 'HO 



(64) 



(65) 



AF^(2) (r; uji, UJ2) 



2/)'2 I , ,2/)/2\ 



in 



\0\0', {u;l0'^+u:l0'^) 



(66) 



13 



Recall that the two Newtonian equations of motion for each inverted harmonic oscillator are given by, 



dr2 



0, Vj = 1, 2. 



Hence, the asymptotic behavior of such macrovariables on manifold M j/jq is given by. 



e, (r) S,e"^^ e M, Vj = 1, 2. 



Substituting 6*1 (t') = Sie'^i^' and 6*2 (r') = 226"^^' into we obtain 



AV" (2) (t; wi, CJ2) 



By direct computation, we find the average of is given by. 



(aK,(2) (r; wi, W2)) 



6 V 



Assuming as a working hypothesis that ^1 = ^2 = 2, we obtain 



ij^g(^i+^2)r' (^S2g2^ir'^2 ^ 2: 



^2„2w2t' ,2 



Finally, substituting (ffTj) in ((64|) . we obtain 



' ^S6a;S^,ifc.i=a;2, 
dr'= J ^sV2S^,if^i»u;2, 



L 18^ ^^2 7 , II OJ2 » UJl. 



2UJT, if UJi — UJ2, 



^M?).^ (r; wi, a;2) oc ■( wir, if wi > ^2, 

W2r, if tj2 > ti^i- 



(67) 



(68) 



(69) 



(70) 



(71) 



(72) 



It is clear that the information-geometrodynamical entropy S ..(2) (r; uji, UJ2) exhibits classical linear behavior in the 

'HO 

asymptotic limit, with proportionality coefficient fl ^ lji + uj2, 

S ^(2) (r; wi, UJ2) ^T. (73) 

Equation (j73p expresses the asymptotic linear growth of our information geometrodynamical entropy for the IHO 
system considered. This result (for n = 2) extends the result of Zurek-Paz (I46p . This result, together with the 
authors previous works [lol . [Tlj lend substantial support for the IG approach advocated in the present article. 

VIII. FINAL REMARKS 



A Gaussian ED statistical model has been constructed on a GA^-dimensional statistical manifold Mg. The macro- 
coordinates on the manifold are represented by the expectation values of microvariables associated with Gaussian 
distributions. The geometric structure of AAg was studied in detail. It was shown that M.s is a curved manifold of 
constant negative Ricci curvature — 3A^ . The geodesies of the ED model are hyperbolic curves on A^^. A study 
of the stability of geodesies on Ms was presented. The notion of statistical volume elements was introduced to 
investigate the asymptotic behavior of a one-parameter family of neighboring volumes ^Vm^ (A) = {Vms {''"i A)}agk+- 
An information-geometric analog of the Zurek-Paz chaos criterion was suggested. It was shown that the behavior of 
geodesies is characterized by exponential instability that leads to chaotic scenarios on the curved statistical manifold. 
These conclusions are supported by a study based on the geodesic deviation equations and on the asymptotic behavior 
of the Jacobi vector field intensity J Ma O'^ M-s- A Lyapunov exponent analog similar to that appearing in the 
Riemannian geometric approach to chaos was suggested as an indicator of chaoticity. On the basis of our analysis 
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a relationship among an entropy-like quantity, chaoticity and curvature is proposed, suggesting to interpret the 
statistical curvature as a measure of the entropic dynamical chaoticity. 

The results obtained in this work are significant, in our opinion, since a rigorous relation among curvature, Lyapunov 
exponents and Kolmogorov-Sinai entropy is still under investigation fs*]. In addition, there does not exist a well defined 
unifying characterization of chaos in classical and quantum physics [5J due to fundamental differences between the two 
theories. In addition, the role of curvature in statistical inference is even less understood. The meaning of statistical 
curvature for a one-parameter model in inference theory was introduced in [s^. Curvature served as an important 
tool in the asymptotic theory of statistical estimation. Therefore the implications of this work is twofold. Firstly, 
it helps understanding possible future use of the statistical curvature in modelling real processes by relating it to 
conventionally accepted quantities such as entropy and chaos. On the other hand, it serves to cast what is already 
known in physics regarding curvature in a new light as a consequence of its proposed link with inference. 

As a simple physical example, we considered the information-geometry M\fjQ associated with a 2-set configuration 
of inverted harmonic oscillators. It was determined that in the limit of a flat frequency spectrum (wi — uj2 — cj), 
the scalar curvature R..{2) (uji, W2) is constantly negative. In the case of distinct frequencies, i.e., uJi 7^ uj2, it 

is possible - for appropriate choices of uji and uj2 - to obtain either negative or positive values of R.^(2) (lui, 102). 

(2) . 

Moreover, it was shown that Mj/jq is an anisotropic manifold since the Weyl projective curvature tensor has a 
non- vanishing component W^i2i2. It was found that the information geometrodynamical entropy of the IHO system 
exhibits asymptotic linear growth. This IHO example is generalized to arbitrary values of n in the Appendix. 

The descriptions of a classical chaotic system of arbitrary interacting degrees of freedom, deviations from Gaussianity 
and chaoticity arising from fluctuations of positively curved statistical manifolds are being investigated . The work 
here presented is shown to be useful to investigate chaotic quantum spectra arising, for instance, from the Poisson and 
Wigner-Dyson quantum level spacing distributions [H, [33, [l^l ■ We remark that based on the results obtained from 
the chosen ED models, it is not unreasonable to think that should the correct variables describing the true degrees 
of freedom of a physical system be identified, perhaps deeper insights into the foundations of models of physics and 
reasoning (and their relationship to each other) may be uncovered. 
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X. APPENDIX 



A. The n-set of inverted anisotropic three-dimensional harmonic oscillators 



We now generalize the results obtained in this article for the n-set of IHOs. The information metric on the 3n- 

(3n) 



dimensional statistical manifold M jjjo given by 
where 

3n 



hn)] ■ 5ij {9i,...., 03„) , 



$(01 



73n 



i=i 



The information geometrodynamical entropy S ..(3^) (r; wi,.., wsn) is defined as 

^IHO 



dcf 



'S'^(3") (t; wi,.., W3„) = lim log [AV (3^) (r; wi,.., W3„) 



where the average volume element AV^(3n) is given by 



(74) 



(75) 



(76) 



T 

(AV" (3^) (r; wi,.., a;3„)) = - / AV^(3,.) (r'; Wi,.., a;3„)dr'. 



(77) 
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and the statistical volume element AV,.(ar,) is defined as 



Substituting ([77)1 and ([75]) in ([75]) we obtain the general expression for Sj.{3,i) (r; wi,.-, f-^sn) 



'S'wO.x) (r; wi,.., W3„) =^ lim log \ - 



To evaluate ([7^ we observe AV (3,1) can be written as 



dr' 



(78) 



(79) 



1 1 / 3ji 

— ^ ( ne' 



3n 



Since the rt-Newtonian equations of motions for each IHO are given by 



dr 

the asymptotic behavior of such macrovariables on manifold M^^^q is given by 

0j (r) Sj-e'^^'^, e M, Vj = 1,..., 3?i. 



We therefore obtain 



Ay (3„) (r; wi,...., uj3n) 



Upon averaging ([83| we find 

AV (3„) (r; wi,...., uj3n) 

3n 



'--3;^ nS, .exp(f]r') 
3n 2~ \«=i 



i=i 



3n 



(80) 



(81) 



(82) 



(83) 



dr'. (84) 



where $7 = '^^^i- As a working hypothesis, we assume = = S Vi, j = 1,.., 37i. Furthermore, assume that 

4 = 1 

ft — !■ cxD so that the spectrum of frequencies becomes continuum and, as an additional working hypothesis, assume 
this spectrum is linearly distributed, 



f^cut-oft 



p{oj) ^ uj with j p{uj)duj = 1, ricut-off = C^, <^ G 



(85) 
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Therefore, we obtain 



> , 1 1 fi ^£^n^\ ' exp(|<17r) 

Ay^(3., (r; ci,...., c.3„) « :;-r3;rS^" ^ (86) 



r — ^oo 



^iHo It 3n 2 

Finally, substituting into (ffS]). we obtain the remarkable result 

3n 



S ^(z-n.) (r; wi,.., LJ3„) ""0:°° f7r, f7 = T^Ci^i. (87) 



Equation (|87|) displays the asymptotic, linear information geometrodynamical entropy growth of the generalized n- 
set of inverted harmonic oscillators and extends the result of Zurek-Paz to an arbitrary set of anisotropic inverted 
harmonic oscillators [J]. 
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